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t-h Abstract 
O 

We investigate five-dimensional braneworld setups with broken Lorentz 
invariance continuing the developments of |1|, where a family of static 
.£h self-tuning braneworld solutions was found. We show that several 

known braneworld models can be embedded into this family. Then we 
give a qualitative analysis of spectra of field fluctuations in backgrounds 
with broken Lorentz invariance. We also elaborate on one particular 
model and study spectra of scalar and spinor fields in it. It turns out 
that the spectra we have found possess very peculiar and unexpected 
properties. 



1 Introduction and Overview 



In modern cosmology and particle physics we are faced with several problems like the cos- 
mological constant problem and the hierarchy problem. Some time ago the novel braneworld 
programme was proposed \2-4\ which made the above problems more tractable |5| (see also |6 
for string theoretical analysis). In this treatment our Universe is represented as a four- 
dimensional submanifold (a brane) embedded into a manifold of higher dimension (a bulk). 
The presence of extra dimensions changes cosmo logical evolution |7| and could be made man- 
ifest at forthcoming high energy experiments |8| and in the other measurements |9|. Among 
an enormous number of the large extra dimensional models the Randall-Sundrum models 
(RSI |3| and the RS2 |4|) appeared to be very successful, as they allow one to look at the 
above two problems from a new viewpoint without spoiling four-dimensional physics at ob- 
servable distances and energies. The RS models have initiated the extended and prospective 
programme of large extra dimensional physics; it is not easy to present an exhaustive list of 
publications devoted to physics in extra dimensional setups, one can use various reviews and 
notes for finding necessary literature |10|. 

In the Randall-Sundrum models the five dimensional anti de-Sitter space is used as a bulk 
filled by the negative cosmological constant and one (RS2) or two (RSI) branes represented 
by 3+1 dimensional hypersurfaces embedded into the bulk. One can also consider various 
generalizations of this configuration and study their phenomenological manifestations. One 
of the possible extensions of the RS models consists of an inclusion of an energy density 
(matter) in consideration which is somehow distributed in the bulk and/or on the brane. 
The presence of matter may lead in general to breaking of Lorentz invariance. Different 
aspects of Lorentz invariance violation in models with large extra dimensions were considered 
in [?, 11 1 . In these models Lorentz invariance was broken spontaneously on the brane due to 
the presence of a symmetry-breaking potential. However, one can also consider models with 
broken Lorentz invariance in the bulk and study deviations from the four-dimensional physics 
on the brane which a priori preserves Lorentz invariance (see, e.g., |12,13,1|). It is worth 
noting however that, as it was proved in |1|, the matter located in the bulk which produces 
a Lorentz violated static smooth metric, inevitably violates the null energy conditions. 

In the present work we investigate field fluctuations in the backgrounds found in |1|, 
where the bulk metric explicitly breaks Lorentz invariance and has the following form 

ds 2 = e~ 2a ^dt 2 - e" 2 ^W - dz 2 . (1.1) 

Here (t, x) are four-dimensional coordinates parallel to the brane and z is the coordinate 
along the extra dimension. In 1 a two-parametric family of solutions of the bulk Einstein 
equations has been found in the presence of an ideal anisotropic relativistic fluid in the bulk. 1 
The functions a and b from (1.1) are linear functions of z 

a(z)=£k\z\, b(z) = (k\zl k>0, (1.2) 

with the coefficients £ and ( related to the anisotropy parameters of the ideal fluid. In this 
paper we introduce a space parametrized by £ and ( which describes different metric tensors 
and study how the fluctuation spectra of fields behave at different values of £ and 

^n accordance with the theorem proved in 1 , this fluid violates null energy conditions. 
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First, we perform more generic analysis of the spectra employing the metric (1.1) without 
specifying the functions a and b. One of our main results is that the spectrum pattern is 
governed by an asymptotic behavior of the functions a and b. If a(z) jb{z) — > as z — > oo then 
the spectrum is discrete and all excitation modes are localized in some vicinity of the brane 
located at z = 0. The reason of the localization is that in this case the potential responsible 
for the localization grows away from the brane. In the opposite case, if a(z)/b(z) does not 
vanish at infinity, then at least part of the spectrum becomes continuous and excitation 
modes with energy high enough delocalize from the brane. In this case the potential tends to 
some constant value (or vanishes) far away from the brane. We further embed some known 
models: flat Minkowski space, the RS2 model |4|, Dubovsky model |12|, into our parametric 
solution. 

In addition, we elaborate on a new model with the metric (corresponding to £ = 0, ( = 1 



where goo metric coefficient is unwarped, and g^- coefficients are warped in the same way 
as in the RS model 2 . The model with metric (1.3) is shown to have dramatically different 
properties than the models |4,12| as long as localization of perturbations and speed of particle 
propagation are considered. In the former field perturbations are localized near the brane and 
their spectrum is discrete, while for the latter there may exist only several (quasi-) localized 
modes and the spectrum is continuous. We study spectra of scalar and spinor excitations 
in the background (1.3). In both cases we find that there is a zero (or gapless) mode — a 
state with vanishing energy as three-dimensional momentum tends to zero. From the four- 
dimensional point of view a zero mode corresponds to a massless state and this is the only 
state which presents in the low energy effective theory. It is remarkable that in spite of the 
metric (1.3) is Lorentz invariant on the brane (z = 0) the zero mode dispersion relation 
explicitly breaks Lorentz invariance at low momenta because of the "tail" of the zero mode 
wave function "feels" the breaking of Lorentz invariance in the bulk. 

The paper is organized as follows. In Sec. 2 we make some qualitative remarks about field 
fluctuation spectra in different backgrounds with broken bulk Lorentz invariance. In Sec. 3 
and Sec. 4 we investigate spectra of scalar and spinor fields respectively in these backgrounds. 
Our conclusions and final comments are given in Sec. 5, Detailed derivation of the fermionic 
spectrum is presented in App. [A] 

2 Localization and Derealization of Fluctuations 

In this section we show that the asymptotic behavior of the functions a and b in (1.1) is 
responsible for the pattern of perturbation spectra of fields in the backgrounds in question. In 
the next sections we investigate how different perturbations behave in a particular example. 
From qualitative description one can understand something about those perturbations in a 
generic fashion. 

2 Five-dimensional coordinates x A are labeled by capital Latin indices 0,1,2,3,5; three-dimensional coor- 
dinates x l are labeled by lower Latin indices 1,2,3. 



in (1.2)) 



ds 2 = dt 2 -e- 2k ^dx 2 -dz 2 , 
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Considering some bulk matter field one can think of (quantum) particles propagating in 
the bulk. Let us assume that our model admits semi-classical approximation. If so, then 
particles do not considerably deviate from their classical trajectories and become "local- 
ized" in domains with higher density of these trajectories. Thus investigation of a geodesic 
configuration helps us in understanding of localization properties of field fluctuations. For 
simplicity we consider a massless free field, so the corresponding particles propagate along 
null geodesies. A geodesic is a solution to the following equation 

x A + Ti c x B x c = , (2.1) 

where Tg C are Christoffel symbols for the metric (1.1) and the dots above x A stand for 
derivatives with respect to the affine parameter. Equations for the 0, 1, 2, 3 components can 
be integrated straightforwardly 

x° = ae 2a , x l = (3 l e 2b , i = 1,2,3, (2.2) 

where a and f3 l are the integration constants. Using this fact and light cone equation for the 
massless particle gABi A i B = one finds the particle velocity along fifth direction 

z 2 = a 2 e 2a - ((3f e 2b . (2.3) 

It follows from this equation that if the three-dimensional initial velocity vanishes, $ = 0, 
then the particle motion is always non-finite. From the quantum mechanical point of view 
it means that the particle with zero three-dimensional momentum cannot be localized near 
the brane. The particle motion is also non-finite in the Lorentz invariant case, a(z) = b(z). 
It does not, however, contradict the fact of presence of a zero mode in the RS2 setup — the 
semi-classical description is certainly not applicable for a zero mode. 

Let us consider the case a(z) — > +oo and b(z) — > +oo as z — > oo. This is the most general 
case corresponding to the finite volume extra dimension. In this case if a(z)/b(z) — > as 
z — > oo then the particle motion is finite (j3i ^ 0) which corresponds to localization of the 
quantum particle in some domain near the brane. The opposite case a(z)/b(z) — > oo or 
a(z)/b(z) — > const as z — > oo corresponds to the derealization regime and there exists a 
probability of a detection of the particle far away from the brane. 

We would like to mention again that the considerations used in this section are merely 
qualitative. An advantage of this description is that its conclusions are universal and can be 
applied to any field theory in the bulk. 

3 Scalar Fluctuations 

After qualitative analysis in the previous section let us proceed to more specific treatment and 
investigate fluctuations of a scalar field in the bulk. For simplicity in this section we consider 
massless scalar field. First, we make a qualitative description of how the perturbations 
behave for generic form of the metric (1.1) and discuss some special points in the space 
parameterized by £ and (. In the end we elaborate on the spectrum of the scalar field in the 
background (1.3). 
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Let us consider the massless scalar field with the five-dimensional action 

+00 



S = J dtdx J dz^g AB d A <\)dB<\>- 



We assume that presence of the field in the bulk does not affect the background metric. 
The equation of motion reads 

[-d 2 + (a' + 3b')d z + e 2a ^d 2 - e 2b ^d 2 } = 0, 

and, after four-dimensional Fourier transform 

it takes the following form 

[d 2 z - {a' + 3b')d z + E 2 e 2a ^ - p 2 e 2b ^] = 0, 
where we denoted p = y/p 2 . For convenience we introduce another function 

X(z) = exp {-\a - |6) <f>{z) , 
and the above equation can be rewritten as follows 

a" + 3b" (a' + 36') 2 1 



d 2 + E 2 e 2a(z) _ p 2 e 2b(z) + 

2 4 

We also introduce the new coordinate y by the relations 

z 

d , 4 



x = o 



dz 



VK Z ) 



3 o(«) 



y(z) 



a{p) d P . 



In terms of this new coordinate, the equation (3.1) takes the following form 

d 2 x 



+ (E 2 -V) X = 0, 



where 



"~KS ,+ Ki) ,+ ' VM 



'd 2 a d 2 b 
l_ 3 

2 \dy 2 dy 2 



(3.1) 



(3.2) 



(3.3) 



Let us consider the potential V in more detail. In what follows we assume Z 2 bulk 
reflection symmetry z —>■ — z with respect to the brane. It means in particular that the 
functions a and b in (1.1) depend on the absolute value of z (and, hence, on y) only. The 
potential V is bilinear in the first and linear in the second derivatives, hence in general it has 
a delta-function term emerging from the last term in (3.3). The sign of the last term can be 
either positive or negative depending on the behavior of the functions a and b in the vicinity 
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of the origin. If the sign is positive then we have a delta-well in the potential and a delta-peak 
otherwise. If Lorentz invariance is not broken (a(y) = b(y)) then the momentum-dependent 
term is merely a constant and does not affect the shape of the potential V. However, if this is 
not the case, depending on the difference b — a in the exponent the shape of the potential can 
be drastically modified. If b > a as y increases into the bulk, then the momentum-dependent 
term increases as well and soon will prevail the other terms. 3 The opposite configuration 
with b < a makes (far enough from the brane) the momentum-dependent term negligible 
compared to the others. Thus we see that the generic potential (3.3) can be classified by the 
following properties 

• Behavior at infinity which is controlled by the Lorentz invariance violation (the sign 
of b — a). The potential V increases/decreases as y — > oo. 

• Sign of the delta-function term. The potential may have either delta- well or delta-peak 
at the origin depending on this sign. 

This rough classification enables us to describe the behavior of fluctuations in each case. 

• If the momentum-dependent term increases as y — > oo the potential has a box-type 
shape and the excitation spectrum is discrete. The potential may have local minima 
and maxima, but qualitative pattern of the spectrum is determined by its behavior at 
infinity. On the contrary, if the potential decays at infinity, then we have continuous 
spectrum of plane waves propagating along y-direction. Some combination of these 
two scenarios is possible if V — > = const as z — > oo. Then those modes with E 2 < 
Voo belong to discrete spectrum and modes with E 2 > Voo contribute to continuous 
spectrum. 

• The sign of the delta-function term affects zero mode existence. In a delta-well there 
might be a zero-modej | and its existence is very unlikely in a configuration with a 
delta-peak. 

Unfortunately rigorous analysis of (3.2) with arbitrary functions a and b is troublesome 
as the corresponding equations cannot be solved analytically. Nevertheless, the qualitative 
behavior of perturbations is not changed if we consider metrics, different in general, but 
having similar form in the vicinity of the origin and at infinity. In Sec. 3.2 we choose the 
functions a and b to be of the form (1.2) which enables us to obtain exact solutions for 
certain configurations. 

3.1 Parameter Space 

In 1 1 1 the static solution for the ansatz (1.1) has been found in the presence of an anisotropic 
perfect fluid in the bulk. The metric has the following form 

ds 2 = e~ 2 ^dt 2 - e- 2 ^lcZx 2 - dz 2 . (3.4) 

3 We do not pay much attention here on rigorous analysis of conditions how the functions a and b affect 
the potential, leaving it to the next subsection where an explicit example is elaborated. 

4 Thcrc might be also no zero mode, but one cannot understand whether it exists or not without detailed 
analysis. 
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Obviously, the chosen parametrization with £, ( and k is overcomplete, 5 one of the parameters 
can be scaled away. Nevertheless, in what follows we keep all three parameters unconstrained 
in order to simplify further analysis. 

One can depict the (£, Q parameter space as shown in Fig. [0 One has flat Minkowski 



< 

c 

III 


1 / 

/ 11 

/ B 




° 

IV 5 



Figure 1: The parameter space. We discuss our theory in the following subsets of this square: 
Minkowski (the origin O), AdS§ (the diagonal), type I (upper triangle) and type II (lower triangle). 

metric at the origin O, AdS$ metric on the main diagonal, the models with £ < ( (type II 
models) in the lower triangle and models with £ > ( (type I models) in the upper triangle. 
For instance, the RS2 model lives on the diagonal, the model proposed in |12| lives at the 
point B. 

For the metric (3.4) the potential in (3.2), (3.3) takes the following form 

V(y) = P \l+ e*|y|) 2K/C_1) + ( f + - (C + 3C)**(l/) • (3.5) 

4(1 + £%|) 

We see that the Lorentz violation contributes the first term of the potential V, which is 
momentum- dependent. The ratio £/£ determines the behavior of the potential at infinity. 
In our considerations we restrict ourselves to the case with "finite volume" of the fifth 



3 For instance, one could impose the constraint £ 2 + Q 2 — 1 and let k play a role of curvature scale. 
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Figure 2: The effective potential for different values of The curves are ordered by steepness: 
for the steepest one the ratio £/£ is 0.3, then the values 0.5, 0.7, 0.8, 1 follow. The brane is located 
at y = and p = k = 1 . 

dimension 6 with the following condition satisfied 



which in our parametrization means that £ + 3( > 0. The corresponding domain in Fig. 1 
has four regions I, II, III, and IV so let us briefly discuss the form of the potential for all 
these regions. First, we assume that both £ and ( are nonnegative (regions I and II). If 
£/£ > 1 then the first term in (3.5) increases at large y thereby producing infinite walls to 
the potential V. Steepness of these walls is controlled by the same ratio The case with 
£ = £ corresponds to the Lorentz invariant background of the RS2 model with AdS^ as a 
bulk space. The case with £/£ < 1 does not differ from the RS2 model qualitatively — the 
potential has the same kind of a shape. In the region III the potential (3.5) has poles at 
\y\ = — (£A;) _1 but they are not physical, our coordinate system is not appropriate for this 
region. However, upon choosing a proper coordinates one can proceed with this region as 
well and observe similar behavior as in the region I. As far as we can conclude from (3.5), 
models located in the region IV do not qualitatively differ from models of type II, indeed, 
the potential has similar asymptotic behavior at infinity. At Fig. 2 we depicted the potential 
(3.5) at different positive values of the ratio 

Employing the considerations from the beginning of this section we can speculate on 
character of the perturbation spectra of type I and type II models. A full exact solution 

6 If the fifth dimension has an infinite volume then there is no in general well defined four-dimensional 
effective theory - the model appears to be five-dimensional without any localization properties on the brane. 
See some considerations in this directions e.g. at 14 . 
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Figure 3: Shape of the model A potential. We put p = k = 1 here. 



for generic value of £/£ is however not possible, nevertheless approximate analysis leads to 
desired understanding of the spectra of these models. Analogously with the considerations 
from the previous subsection we can straightforwardly state that the models of type I have a 
box-type potential and hence have discrete spectrum of the scalar perturbations, whereas the 
models of type II have decaying potential at infinity and plane wave excitations propagating 
along y-axis. The excitation spectrum of type II models is therefore continuous in general. 
Due to the delta-well it may contain a zero mode as well. Properties of type II world were 
discussed in |12| (point B). Now we shall investigate the scalar fluctuations in the completely 
different background A with almost opposite properties. 

3.2 The Model A. Generic Solution and Static Potential 

In order to proceed with the model A corresponding to the point (0, 1) in Fig. 1 one considers 
the Schrodinger equation in the form (3.1). The potential reads 



spectrum to be discrete. Schematically the potential (3.6) is shown in Fig. 3. First, we solve 
(3.1) with the potential (3.6) for z ^ 



then we impose the boundary conditions at the origin to take into account the delta-function 
term 



V(z) = p 2 e 2kVA + -k 2 - 3k5{z) . 



(3-6) 



We see that the potential tends to infinity as \z 



oo, thus we expect the perturbation 





(3.7) 




(3.8) 
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where v 



| — f§ and N is a normalization constant which is determined by the condition 



+00 



dz\ X (z)\ 2 = l. 



(3.9) 



The matching condition (3.7) at z = yields 

P*W(t) 3 



fc K„ (f ) 2 



(3.10) 



This equation relates the energy and momentum of modes, therefore it represents the dis- 
persion relation. We consider it below. 

Green Function. Let us now find the brane-to-bulk Green function for the scalar field (p. 
In the momentum space it satisfies the following equation 7 

-8 2 - E 2 + p 2 e 2k ^ + ^k 2 - 3k5(z) A p (E,p,z) = 5(z). (3.11) 

Using (3.8) one obtains the decaying solution to this equation outside of the brane 

A p (E,p,z) = C(E,p)K u ge^l) , 

where C(E,p) is yet to be determined constant. Due to the presence of the delta-functions 
in (3.11) one imposes the following constraint 



[d z A p {E,p,0)\ + 3kA p (E,p,0) 



-1. 



This equation determs the constant C(E,p). Now we are ready to find the propagator. Let 
us put z = in order to extract the brane-to-brane propagator 



AJE,p,0) 



k 



IW (j) 

k MI) 



p 



V 



We see that the condition for the brane-to-brane propagator to have a pole exactly coincides 
with (3.10). Later we show that (3.10) has only real valued solutions in terms of E, hence 
all perturbation modes in our model are stable. In order to find the static brane-to-brane 
Green function we put E — and obtain it in the momentum space 



G p (0) = A p (0,p,0) 



1 

2k 



p K|(g) 



-1 -1 



- 3 



k 1 
2^2 + 2p 



(3.12) 



7 Due to SO(3) rotational brane symmetry the Green function depends only on absolute value of momen- 
tum p. 
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Inverted Fourier transform leads us to the Green function in the coordinate representation 
which represents the static potential 



+00 

<f'l> ■ 4tt f (k \ ,>'/"■-<■-'/"■ 



G ( r ) = / 7^Gp(^ pr = TT^r dp[- + l 



(2?r) 3 pv 1 (2vr) 3 r J \p ) 2i 



+00 +00 



fc / dp — ^ - 2i / dpe ipr = 1 + — . (3.13) 



2ir 2 r 1 J p J J Anr \ nkr 







The obtained potential describes the Coulomb (Newton) law with small correction at large 
distances compared with the curvature scale 1/k. At small distances the interaction becomes 
five-dimensional. Note that unlike the RS2 model |4| the correction to the potential has a 
different power: r~ 2 instead of r~ 3 for the RS2 model. 



3.3 Spectrum of the model A 

Let us now discuss the spectrum of the scalar sector in the model A. We proceed with 
expanding (3.10) at small and large momenta and derive dispersion relations for different 
excitation modes. The spectrum of the scalar field including the zero mode and higher 
modes is shown in Fig. 4, We observe that there exists zero mode and higher excitations. 
The behavior of the higher excitations at small momenta is of particular interest, as all 
curves collapse to one point E = \k as p — ► 0. 



Zero mode at small momenta. From (3.10) one observes that there exists a solution 
with E(p) —>■ as p —>■ which corresponds to the zero mode. The corresponding dispersion 
relation can be found if one performs an expansion in the vicinity of the point E — 0, p — 0, 
and after straightforward calculation one obtains 

E 2 = 3 p 2 (l - I + 0{p 2 )) . (3.14) 

We see that the zero mode dispersion relation is not Lorentz invariant. However, the low 
energy theory is approximately Lorentz invariant in the sense that the dispersion relation 
E 2 pa 3p 2 has the form of the ordinary relativistic formula with rescaled speed of light. From 
(3.14) one finds the speed of propagation of scalar particles (group velocity dE/dp) which 
is equal to \/3 and exceeds a speed of tightly bounded massless particles on the brane being 
equal to unity. It also exceeds the speed of scalar particles considered in |12| which is always 
less than 1. Therefore, one can expect that while moving in the parameter space in Fig. 1 
from point B to point A throughout the upper right region, the speed of scalar particles will 
increase from the values found in |12| to y/3 derived here. 

Let us calculate the zero mode contribution to the static brane Green function (3.12) 

«f(0) = «. (3,5) 
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In order to find it one should know the zero mode normalization constant Nq which can be 
determined at low momenta in the following way. One notes that dK u (x) jdv = K u (x) (log x+ 
(9(1)) (see, e.g. |15|), so the zero mode wave function can be represented as 



N K 



k 



l + O p^logfe 



After that the normalization constant can be calculated from (3.9) 



(3.16) 



m = k- 



71 



o 



¥ log k 



Given Nq and substituting the zero mode dispersion relation (3.14) together with (3.16) into 
(3.15) one finds the zero mode contribution to the static Green function 



G(°)(0) = A + ^ + O(logp) . 



(3.17) 



We see that at low momenta the Green function (3.13) is saturated by the zero mode. The 
corrections 0(logp) emerging in (3.17) should be canceled out by contributions from higher 
modes. 



Higher modes at small momenta. As was mentioned above the mode equation (3.1 ) can 
be written as a Schrodinger equation with a potential dependent on a particle momentum. 
We can apply WKB method here for understanding the spectrum of higher excitations 8 . The 
quantization law reads 9 

j ^El-V{z) dz = 27m, (3.18) 

where the potential V(z) is given in (3.6). Because of Z 2 reflection symmetry of the potential 
we have 

2-2 J y/E* - V{z)dz = 2vrn , (3.19) 
o 

where the classical turning points have the following coordinates for each n 

2k p 2 

Calculating integral in the l.h.s. of (3.19) one has 



El - \k 2 - p 2 + yj El - \k 2 Arctanhy i - — - gfc2 = —n 



p 2 ixk 

V 



8 In order to find the dispersion relation for the higher modes one can solve (3.10) at the low momenta 
directly. However, we do not perform this analysis here. For the fermionic sector, where the spectrum of 
higher modes will appear to be very much similar to the scalar one discussed here, we shall expand an analog 
of (3.10) directly. 

s In our calculations we ignore ^-function contribution to the potential V as it brings only small correction 
to the integral in the l.h.s of the WKB formula. Therefore we can also neglect | in the r.h.s. as it is irrelevant 
for our precision. 



13 



Expansion for < p <C k yields 
or 

^ 9 9, 9 n 2 k 2 n 2 

E »=4* 2 + ^' (3 ' 2 °) 

One can show that the WKB method used here is applicable provided that p/k < e~ n . 

The obtained result (3.20) becomes more transparent if one notes that at small momenta 
the potential (3.6) can be approximated by the box-type potential 



1 p 

if \z\ < — — log - 
k k 

1 p 

oo if \z > — — log — 



Then the dispersion relation (3.20) immediately follows. Note also, that if p = the spec- 
trum becomes continuous what corresponds to a particle which can propagate to infinity in 
complete agreement with the results of Sec. 2, 

Large momenta. At large momenta p ^> k one can neglect the factor of | in the r.h.s. of 
(3.10) and (3.10) becomes equivalent to the following condition 

k< m = o. 



V 



This condition is considered in App. A in the context of fermionic spectrum at large momenta 
(see (A.8)). The result is 

E = p + 0{p 1 z). 

Thus at large momenta, as expected, the theory becomes Lorentz invariant. 



4 Fermionic Fluctuations 

In this section we discuss the spectrum of fermionic fluctuations in braneworlds with broken 
Lorentz invariance. One can perform a qualitative analysis for generic (£, Q configuration 
(3.4) and reduce the appropriate five dimensional Dirac equations to a Schrodinger-type 
equations. Then, analogously to the scalar sector discussed at Sec. 3, one can establish that 
the asymptotic behavior of the metric at infinity is responsible for the actual pattern of the 
spectrum. For instance, if £ > ( the spectrum is continuous. In this section we refrain 
from considering qualitative picture with generic ansatz (3.4) but focus our attention on the 
model A where exact solution is obtained and analyzed. 

We consider the four-component spinor field ^ on the five-dimensional manifold with 
metric (1.3) 

+oo 

S= dtdx / dzy/g^f-® , (4.1) 
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where is the five- dimensional mass of the fermion which may however depend on z. In 
this section we study two cases: 

i) m^p = m ■ sign(z). It corresponds to a usual (Lorentz invariant) way of localization of 
fermions in extra dimensional setups (see, e.g. |16| and for more recent treatments |17|). It 
is worth noting that in this case the action (4.1) is invariant under Z 2 reflection symmetry 
z — > —z if the fermion field transforms as 

*(-*) = Tl^(z), (4.2) 

where phase =F can be only determined by fermion interactions. 

ii) = m. In this case we deal with a usual massive fermion and the action does 
not possess the Z 2 symmetry. It does not play a role, however, since this case has only an 
illustrative character and we are interested in the limit m — > which can be easily obtained 
in this case and which restores Z 2 symmetry. 

Generic Solution. The Dirac equation obtained from (4.1) reads 

iT A V A ^{x,z) -m^(x,z) = 0, (4.3) 
where Va = <9a + oja is the covariant derivative with the spin-connection 

u A = (0, -|A;757ie" fc|21 , 0) sign(z) . 
The five matrices T A (A = 0, 1, 2, 3, 5) obeying the Clifford algebra relations 

{T A ,T B }=2g AB , 
are related to the 7-matrices (/i, u = 0,1, 2, 3) 

{^, Y} = 2iT , 7 5 = 75 = ^7°7 1 7 2 7 3 , (4.4) 

as follows 

r° = 7 , p = e fc|2| 7\ r 5 = -r 5 = -i 7 5 . 

In order to simplify further calculations it is convenient to work with the rescaled field 

# = et fc|2| V, (4.5) 
and apply the four-dimensional Fourier transform to the field ip 

V>(t, x, z) = J dEdp e lEt - l ^(E, p, z) + h.c. . 

In this relation and in what follows we assume that E > 0. Then due to (4.5) the normal- 
ization condition for the field ip(E,p,z) has the canonical form 

+00 +00 

J dV#*7o*= J dz^(E,p,z)^(E,p,z) = l. (4.6) 
—00 —00 
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In terms of the new field the Dirac equation (4.3) takes the form 

(£ 7 ° - e fcN 7 3 P - l 5 d z ) i> + m^ip = . (4.7) 

where we set p = (0, 0,p) for simplicity. In order to solve this equation we use the following 
decomposition of the spinor 

tj) = ^<j> a ji{E,p,z)U 0h p, (4.8) 

where a and (3 independently take values ±, (fc a ^ are four unknown functions, and U a% p are 
four constant spinors with the following properties 

2S 3 U ai p = i7i7 2 t/a,/3 = otU a fi 
loU a ,f3 = a(3U a ^ 

lsU a , p = -U a> - P . (4.9) 

Let us say a few words concerning the above relations. Because of S3 (the third component 
of the spin) commutes with (4.7) it is easy to use the basis composed of the eigenvectors of 
S3 — this is what the first property tells us. Since 70 commutes with S3 one can choose 
the basis in which the eigenvectors of S3 are eigenvectors of 70 as well — this is what the 
second property tells us. The first two properties correspond to fixing of the orthogonal 
basis up to four unknown constants. Since 73 commutes with S3 and anticommutes with 70 
one can check that 73 ~ U a _p — the third property fixes the constant in this relation 
and, thereby, fixes the basis up to two constants which can be found (up to phases) from the 
normalization conditions 

Ua'fi'Ua^ = 25 a ' a 5i3'i3 . 

The last property follows directly from the first three properties and from the definition of 
75 (4.4). 

Substituting the expansion (4.8) into (4.7) one obtains the following equation 

(d z - /W M )0«,/3 = (a(3E - m^.p . (4.10) 



Acting on both sides of this equation by the operator d z + /3pe kz and using (4.10) with 
(3 — > —(3 one ends up with the second order equation 

(d 2 z +E 2 -m 2 - p 2 e 2fc|z| - (3pk sign(z)e fc|21 ) (j) a ^ = -2m5{z)<j) a ,_p , (4.11) 

where we assume that = m ■ sign(z). It is easy to check that the normalizable solution 
of (4.11) at z > 0, p > has the form 10 



(f e* 1 * 1 ) , (4.12) 
where 

&(z) = [K, + i(x) - PK^iix)] . (4.13) 



■'No summation over repetitive indexes a and /3 is assumed. 
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In what follows we use the superscripts ">" for the functions at z > and "<" for the 
functions at z < 0, C^'^ are constants, and v = k~ x \Jm 2 — E 2 . It is worth noting that the 
Dirac operator in the problem under consideration is Hermitian. Thus the energy E has to be 
real (in fact it follows directly from the dispersion relation). It means in particular that v has 
to be either real [E < m) or pure imaginary {E > m). In both cases the sign of v does not 
play any role: it can be absorbed into the definitions of C a ^ because of £,p(—v) = —P£,p{v). 
So we assume that either real (E < m) or imaginary (E > m) part of v is positive. 

One can see from (4.11) that under the reflection transformation z —>■ —z only the last 
term in the l.h.s. changes the sign. Thus the solution at z < can be represented as 



It follows from (4.11) that the solutions for p < can be obtained from (4.12) and (4.13) by 
interchanging (3 — > —(3. The eight unknown constants C^'p have to be determined from the 
normalization condition (4.6), from the continuity conditions at the origin 

^(0) = ^(0), (4.14) 
and from the matching conditions for the derivatives which follow from (4.11) 

d z <f>l p {0) - d z <p<p(0) = -2m0 a ,_^(O) . (4.15) 

At fixed a, (4.14) and (4.15) form a set of four linear homogeneous equations for the four 
unknown constants C^'p- Hence in order this set to be self-consistent its determinant should 
be equal to zero. Calculating the determinant we end up with the equation 



/ (9 \ 2 Am 2 

t 2 (^(e+(t)e-(t))J =-^(£ + e-) 2 , (4.i6) 

where t = p/k. This equation represents the dispersion relation E = E(p) for the fermion 
field in question. 

The equations (4.15) and (4.16) can be considerably simplified if one notes that the 
functions ^ satisfy the relations 

—, = — + PQ3 ■ 

ax x 

Using these identities one obtains from (4.14) and (4.15) the relations for the constants C^'g 

Cl- = ip%C>^, (4.17) 

where 7 = ±1 corresponds to the two possible signs of the square root of (4.16). Making use 
of (4.9) one verifies that 7 corresponds to the =F signs in (4.2). Thus 7 remains ambiguous 
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and can be only fixed by the fermion interactions. As it follows from (4.17), at fixed a we 
have only one unknown constant, say 



C a = Ut)C> + = lUt)C>_ , (4.18) 

which is fixed by the normalization condition (4.6). It means that at fixed a and positive 
energy we have only one degree of freedom — a particle with the polarization a. 

It is worth noting here that v cannot vanish at any value of energy. Indeed, if v = 
then £ + = 0. It means that C a _ = (see (4.17)), and, so, both <j) ay ± = 0, i.e. there is no 
normalized solution in this case. 

Let us now substitute (4.17) into (4.12) and afterwards into (4.10). Then one finds that 

LM = , (4.19) 

where we have introduced the dimensionless energy e = E/k and mass \i = m/k. In fact this 
equation does not depend on (3: if one changes (3 — *> —(3 then by making use of the definition 
of v in the r.h.s of (4.19) one obtains the same equation. Thus (4.19) is equivalent to the 
equation 

£_ (t) u — ae / \ 

(420) 

This relation is nothing but the dispersion relation E = E(p). One can check that (4.20) 
follows directly from (4.16). However, from (4.16) one can obtain another branch of solutions 
which corresponds to the substitution a — > —a in (4.20). This branch does not satisfy (4.10) 
and should be rejected. 

Let us now consider the case = m. In this case the term with the delta-function in 
the r.h.s. of (4.11) is absent. Equations (4.12)— (4.14) still hold, in (4.15), (4.16) the r.h.s. 
vanishes. Thus the equation yielding the dispersion relation in this case (an analog of (4.20)) 
follows from (4.16) 

e-(*) 



where 7 = ±1. An analog of (4.17) reads as follows 

y< _ &(*) n > 



ij, (4.21) 



The second equation in (4.22) follows from (4.10). This equation in fact does not depend 
on f3 — interchanging (3 —(3 leads to the same equation if one takes into account the 
definition of v. Thus we have no other constraints on the constants G a $ and at fixed a we 
have only one degree of freedom. 

Let us now discuss the spectrum of the fields <p a ^ which follows from (4.20) and (4.21). 
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Spectrum of the model A 



Here we give a summary of the fermionic spectrum of the model A at low and high momenta 
and at different values of the bulk mass. The detailed calculations are given in App. A, The 
spectrum pattern at low momenta is shown in Fig. 5, We begin with the case = m-sign(z). 



Zero mode at low momenta. A zero or gapless mode is a state with energy E(p) — > 
as momentum p — > 0. From the four-dimensional point of view a zero mode corresponds 
to a massless state. As we will see the zero mode is present in the fermionic spectrum of 
the model both in the case m = and m ^ 0. However, the dependence of the energy 
on the momentum p <C m, k is different for different values of the bulk mass m. Here we 
consider only the case of nonzero mass. As is shown in App. A, the dispersion relation has 
the following form 

From this formula the following three regimes can be derived 
• At m > k/2 the dispersion relation is 

f 2m \ /t ts 

M^r- (424) 

We see that in the limit m 3> k the dispersion relation becomes Lorentz invariant as 
it should be as this limit corresponds to the flat case. 



At m = k/2 we have 



pe lE 



where 7# ~ 0.577 is Euler constant. 
At < m < k/2 we have 



2m 



■p ( 1 TO \ 

E ~ 2m-)| k 4 ( 4 J " ■ (4-26) 



r(| + f) \2k 

In all these regimes it turns out that a = — 1, and 7 = 1. Thus the zero mode describes a 
massless polarized particle with helicity — 1. 

Higher modes at low momenta. We start with the case = m ■ sign(z). Similarly to 
the scalar spectrum discussed in Sec. 3 one can study the spectrum of higher fermionic modes 
employing the WKB approximation. Here we, however, present the dispersion relations for 
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Figure 5: Fermionic spectra e = e{t) for \x = |, 7 = 1 and = m • sign(z) (upper plot), and 
H = (lower plot). Straight line starting from the origin of the upper plot shows the zero mode 
dispersion relation. Other curves originating from the point (0, |) show dispersion laws for higher 
modes at different n = 1, 2, . . . . At the lower plot all curves start from the origin. 
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the higher modes obtained directly from (4.20). Solving this equation at p <C m, k one finds 
(the detailed calculations are given in App. A) 



E„ 



imk 



m \ 1 + tt\ 



2m log — 



1 + 



2k 



irnk 



2m log f 



(4.27) 



where ^f(z) is Euler digamma function and n e N is the mode number. The helicity depends 
on the mode number as a = (— 1)™ +1 . The latter approximate equality in (4.27) is obtained 
by the WKB approximation and is similar to (3.20) for the scalar case. It is worth noting 
that the first equality in (4.27) is valid only for n such that the term in the parenthesis is less 
than unity. At n large enough the character of the spectrum does not change: it remains to 
be discrete and E — > m as p — > 0. 

It follows from (4.27) that at fixed level number n there are two states with different 

Moreover, these two massive states are polarized and 
So the spectrum of the model for the positive energies 



energies E{pf = 1) and E(j = —1] 
have the same helicity a = (— l) n+1 , 
is completely non-degenerate. 

Now we consider the case with the constant bulk mass ■ 
p <^ m,k. In this case we obtain the following dispersion relation 



m at small momentum 



E„ 



\ 



m 2 + 



kn(2n + 1) 



4* (I) 



41og^ 



(4.28) 



2k 



where n is a nonnegative integer number. We also have 7 = (— l) n+1 . It is worth noting that 
7 rather then a is fixed in this case and E n does not depend on a at all. Thus for fixed n 
the spectrum is doubly degenerate - there are two modes with helicities a — ±1. One also 
sees that there is a smooth limit m — > in this case in which all modes become gapless. 



Modes at high momenta. At momenta large enough, p ^> m,k, the dispersion relations 
for all modes in the both cases m^ = m ■ sign(z) and m^ = m become 

E = p + 0(pl), (4.29) 

and the correction in (4.29) is positive. Thus at large momenta the Lorentz invariant dis- 
persion relation is restored. 



Chiral properties of solutions. It is a standard lore that in the Lorentz invariant case 
for m^p = m ■ sign(z) a zero mode appears in the spectrum and this zero mode is chiral from 
the four-dimensional point of view. Let us consider what happens in our case. Using (4.9) 
one finds ^ 



n It worth to note that 7 (sign in (4.2)) plays a role of the parameter of the model. Thus for the fixed 7 
there is only one state. 



21 



From this equation by means of (4.12) and (4.18) one has 



The corresponding expression at z < can be obtained from (4.30) by interchanging <-> £ + . 
Thus both left-handed and right-handed components do not vanish as functions of z. In 
particular, the zero mode contains both left-handed and right-handed spinors. However, for 
the zero mode 7 = 1 and the right-handed spinor vanishes on the brane z — 0. So, in this 
respect, one can say that only left-handed spinor is localized on the brane. This statement 
is true for higher modes as well: if 7 = 1 only left-handed spinors are localized on the brane, 
whereas, right-handed spinors are localized for 7 = — 1. In both cases depending on a the 
modes have different helicities (for the zero mode a = — 1 and helicity is —1). In Fig. 6 some 
mode profiles from (4.30) are shown. 




n = 2 n = 3 




Figure 6: Mode profiles at n = 0,1,2,3 for m = 2k, p = 0.01. Red solid curves correspond to 
left-handed and dashed blue curves correspond to right-handed spinors. On the upper left plot 
(n=0) the right-handed component is strongly suppressed. 
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In the case = m an analog of (4.30) is 

^{^^^)%u M ( , 31 , 

so the chiral components do not vanish even on the brane. However, in the massless limit 
the expression in the parentheses in (4.31) tends to 



Thus depending on 7a either left-handed (7a = 1) or right-handed (7a = —1) modes 
are localized on the brane (again in the sense that the corresponding spinors vanish on the 
brane). In particular, for the lowest modes (n — 0) 7 = — 1 the localized spinor is left-handed 
if helicity is a = —1 and right-handed in the opposite case. Making an insight into (4.31) we 
see that even for the constant mass localization of fermionic perturbations is possible. The 
reason for that is based on the special geometry we work in and on negative curvature of the 
space I IS I . 



Propagation velocity and relation to scalars. To conclude this section let us note 
that from eqs. (4.24) — (4.29) one sees that group velocity dE/dp in all regimes is greater 
than unity. Moreover, it follows from eqs. (4.25), (4.26) and (4.27), the group velocity tends 
to infinity as p — > 0. Thus in order to have equal velocities at low momenta for the scalar 
(c s = a/3) (3.14) and the fermionic zero mode one needs to involve some sort of fine tuning. 
Namely, one assumes that m > k/2 (only in this case the group velocity is finite), uses (4.24), 
and solves the equation 

1 2m 



2m — k 

The solution is ? ~ 1.18 > 1/2 which confirms that our assumption is valid. 



5 Conclusions 

In this work we investigated perturbation spectra of the scalar and the spin-| fermion fields 
in the extra dimensional setups with broken Lorentz invariance in the bulk solutions found 
in 1 1 1 . In the beginning we qualitatively discussed the relation between asymptotic behavior 
of the bulk metric tensor at infinity in an extra dimension and localization of higher (non- 
zero) modes by analyzing geodesic equation in the given metric. Then for the massless 
scalar field in the bulk we performed the spectrum analysis by reducing the corresponding 
mode equation to the Schrodinger equation with a potential dependent on the metric and 
particle momentum. Existence of a zero mode is given rise by the presence of a delta- 
well in a potential. We have constructed the parameter space of the bulk metrics used 
in J 1 1 with explicitly broken Lorentz invariance in the bulk and discussed how spectra of 
excitations behave in different domains of this space. It has been shown that if the goo 
metric coefficient increases faster then ^ as the fifth coordinate tends to infinity, then 
higher modes of field fluctuations are not localized and the spectrum is continuous; if the 
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opposite case is realized then a spectrum is discrete and all modes are localized on the brane. 
Localization properties may, however, depend on three-dimensional momentum. In the end 
of Sec. 3 we have elaborated on the model with goo = 1 and gij = e -2 ^ 2 '^ as contrasted to 
the model |12|, in which g Q0 is warped while g^ is not. 

The static potential for the massless scalar field was derived. Unlike the RS2 case our 
potential acquires different power correction at small distances which could make the two 
models experimentally distinguishable. 

The dispersion relation for the zero mode was obtained and appeared to break Lorentz 
invariance explicitly. The speed of scalar massless particles propagating along the brane 
(which is governed by the zero mode dispersion relation) was found to be larger then the 
speed of tightly bound massless particle on the brane. Recall that in |12| the speed was 
smaller and one can expect, that for intermediate values of Lorentz violating parameters a 
speed of scalar particles will be bounded by our result from above and by the result |12| from 
below. However, for completeness, other domains of the parameter space Fig. 1 have to be 
analyzed. 

In Sec. 4 spinor fluctuations were investigated on the example of the model A. We have 
explicitly found the solution of the Dirac equation with the mass term = m ■ sign(z) 
and derived the corresponding dispersion relations. We found fermionic spectra for different 
values of the ratio m/k. Generically these spectra appeared to be almost the same as in 
the scalar case except for the zero mode dispersion relation. For m > the energy E is 
proportional to the momentum p with the coefficient dependent on this ratio; for m = ^k we 
had E ~ plogp; and m/fc-dependent power law behavior for m < \k. Moreover, it turns out 
that all modes including the zero mode are polarized but nonchiral 12 unlike the usual Lorentz 
invariant case where an appeared zero mode is always chiral. Thus the energy spectrum is 
non-degenerate. From the four-dimensional point of view the nonchiral character of the zero 
mode means that it is not described by the Weyl equation — there is a mixing between 
left-handed and right-handed spinors (mass-like term) which depends on three-dimensional 
momentum and disappears in the limit of small momenta. Apart from that, chirality is 
precisely restored on the brane since either left or right-handed components of spinor wave 
functions vanish on the brane (the same is true for the higher modes as well). 

We have also investigated the case = m. It has an illustrative character due to the 
broken Z2 bulk symmetry. One can take the limit m — > and Z2 symmetry is restored. 
Unlike all the other known braneworld models, a zero mode does exist even for m = which 
can be considered as a consequence of our special geometry and negative curvature of the 
space |18|. Namely, all excitation modes can be treated as zero modes since for all modes 
energy vanishes as p — > 0, unlike the scalar massless case where there is a gap 3k/ 2 between 
zero and higher modes. However, the modes are not polarized and the energy levels are 
doubly degenerate. Moreover, all modes are nonchiral as in the m$ = m ■ sign(z) case. 

Both in the scalar and fermion cases we found that the spectra of higher modes have 
very peculiar properties — at small momenta energies of all modes tend to one and the 
same value independent of the mode number. From the four-dimensional point of view it 
means that all excitations have the same mass. However, on the other hand a particle with 
precisely zero momentum becomes delocalized and may freely escape from the brane. It is 

12 The modes correspond to a mixture of left-handed and right-handed spinors with certain coefficients. 
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worth stressing that for a particle somehow localized in the three-dimensional space (wave 
packet) the momentum cannot be exactly equal to zero due to the uncertainty principle. 

To conclude, the considered model possesses a lot of unusual and remarkable properties 
which are interesting from theoretical as well as from phenomenological points of view. 

Note added. Recently we found an interesting paper |19| where the authors deal with 
metric deformations which break Lorentz invariance similar to those considered in our work. 
In 1 20 1 an explicit comparison of the two solutions is made. The work |19| represents a 
prospective programme of gauge/gravity correspondence for theories with broken Lorentz 
invariance. 
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A Appendix 

In this Appendix we give detailed derivations of the dispersion relations (4.24) - (4.29) for 
the fermionic spectrum. For simplicity we introduce the following dimensionless variables 
t = p/k, e = E/k, and fi = m/k. We start from the equation (4.20). 



Small momenta. Let us firstly discuss the regime of small momenta t C 1, One can 
expand the l.h.s. of (4.20) at small t and obtain the following equation 



t 



+ 



r(i/2 



2v- 1 r(l/2 + i/) \2 



2v 



7- 



ae 



v 



t 



T(l/2-u) ft 

2v- 1 T(l/2 + v) \2 



2;/ 



• (a.i; 



Here we kept two terms in the next to the leading order since depending on v these terms 
may compete with each other; the corrections to the equation are of the order of t 2 and t 2u+1 . 
If Re v ^ (which automatically implies /i ^ 0) then the leading order is represented by 
the first terms in the square brackets in (A.I). Thus in order to satisfy the equation one is 
forced to conclude that e — > as t — > and 7 = 1, which means that we deal with the zero 
mode. Aside from that, it turns out that a is negative. This fact follows straightforwardly 
from (4.20). Indeed, the operator (4.11) with the boundary conditions under consideration 
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is hermitian, so, if Re v 7^ then Imz/ = 0. It is well known that at real values of order the 
modified Bessel function K u (t) is positively defined and at fixed argument is an increasing 
function of v. Thus, > £,+{t) > what means that a = — 1 and 7=1. 

In order to find the behavior of the zero mode at t — > one needs to take into account 
subleading terms in (A.l). Expanding (A.l) at small e and t one obtains 

- 2 " -« + 2*4ar (a.2) 



2/x-l r(| + //) V2 

which reconstructs all dispersion relations for zero mode (4.24) - (4.26) at small momenta. 
Indeed, if /i > 1/2 then the first term in (A.2) is dominant, the second term is dominant 
in the opposite case (// < 1/2). In the case /1 — > 1/2 both terms tend to infinity becoming 
comparable. Expanding the r.h.s. of (A.2) at fi = 1/2 one easily verifies that divergent parts 
are cancelled and obtains (4.25). 

Recall that in the above analysis we assumed Re v 7^ 0. Let us now consider the case 
when v is pure imaginary. In this case the leading terms in the brackets in (A.l) are the first 
and third one. Then after some calculations one obtains from (A.l) the following equation 
(the same equation follows directly from (4.16)) 

t\- 2v (\ \\ 4 (t\ 2v (\ - 2 



r U + Z/ ) (^-^-U) T {2~^ (7A* + ^) = 0. (A.3) 



which is equivalent to 



, t ir(i/2 + i/)-r(i/2-i/) 

2v log — h 2i arctan - — - — ; r — — ; + % arctan 

6 2 iT(l/2 + u)+V(l/2-u) 



7— 



Tmi , (A.4) 



where n is an integer. Since the last two terms in the l.h.s. of this equation are finite at all 
values of imaginary v, the only way to satisfy the equation at fixed n is to assume that v 
thends to zero as t — > 0. Note, however, that the expression (A.4) does not depend on a, i.e. 
it holds for both helicities. In order to clarify the dependence on a one should find v directly 
from (A.l) (the advantage of (A.3) is that this equation allows one to establish quite easily 
the fact that v — > as t — > 0, while it is not so obvious from (A.l)). 
Performing small v expansion in (A.l) one obtains 

tanh^(|)-Wog0 = ^y . (A.5) 

This equation tells us that if a = 1 then the argument of tanh should go to + |), I e Z, 
meanwhile if a = — 1 the argument tends to inl. Then the solution of (A.5) can be written 
as follows 

vn(t) = J , (A.6) 

»(»)- I-21og- 

where n G N (recall that v ^ and we assume that lm(v) > 0), and mode helicities depend 
on the level number as a = (— l) n+1 . The above equation in turn gives us the dispersion 
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relation for the heavy modes (4.27). Note that (A. 6) is obtained in the assumption \u\ <C 1. 
Thus for n large enough the momentum t should be small. However, it does not affect the 
conclusion that the spectrum is discrete and E n — > m as t — > 0. 

Let us now consider the case = m, i.e. we start from (4.21). Comparing (4.21) with 
(4.20) one sees that (A.l) still holds if one replaces the factor (fi — ae)/i> in the r.h.s. by the 
factor of i which is equivalent to taking \x = and a = 1 in this (and only in this) factor. 
The equation (A. 3) also holds if one puts 7/i = in the parentheses. So the conclusion that 
v goes to zero as t — * remains intact, i.e. all modes have the same mass /1, in particular in 
the massless case all modes become gapless. In the regime v — ► one can obtain from (A.l) 
the following equation (an analog of (A.5)) 

tanh [v^l - v log = -17 . 

Then 

. . m{2n + 1) . . 

Vn(t) = rr; —r • (A.7) 

Note that 7 depends on the mode number n as 7 = (— l) n+1 . From (A.7) the dispersion 
relation (4.28) immediately follows. 

High momenta. Let us now consider the high momenta regime t>/i,6. In this case it 
is easy to start from (4.11) and neglect the last term compared with the term proportional 
to p 2 . Then the solutions to the equations of motion reads 

e± = K I/ (t), v= vV-^ 2 , 

and the equation yielding the dispersion relation (an analog of (4.16)) is the following 

tK' u (t)=0, (A.8) 

where we neglected the term fiK u which appears in the case m$ = m ■ sign(^). Thus one can 
consider both cases (m^ = m ■ sign(,z) and = m) on an equal footing. 

There are three different possibilities in the regime t — > 00: 1) \v\ <ti t, 2) \v\ > t, 3) 
\u\ — > 00 but \u\ < t. In the first case one can use the asymptotic behavior of the Bessel 
function with large argument and fixed order 

M*) = ^e-*(l + 0(r 1 )) . (A.9) 

Due to the fact that the leading order in (A.9) does not depend on v there is no solution 
of (A.8) in this case. In the second and in the third cases e — > 00 as t — > 00, so v is pure 
imaginary v — i\v\. In the second case one uses the asymptotic expansion for the Bessel 
function |15 



V2 



I v 1 7r 
2 



r {\) sin (V|Arccosh^ - ^\u\ 2 -t 2 + ^ J + O (r 1 ) 
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Substituting this expression into (A. 8) one obtains the equation 



3 



tan f |i/|Arccosh^ - ^\u\ 2 - t 2 + ^ = 2 ^ - 1 ^ . 
The solution of this equation reads 

\ v \ =t + 0{t 1/3 ), (A.10) 

where the correction is positive, so the condition |z/| > t is satisfied. From this equation one 
obtains (4.29). 

In the last case \v\ < t one uses the asymptotic expansion |15 



r(|) + o(r 1 ; 



\v\" — \v\ arcsm 



Substituting this expansion in (A. 8) one again ends up with the solution (A. 10) where the 
correction is positive as well. It contradicts however the assumption \v\ < t. Thus there is 
no solution of (A. 8) in this case. 
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